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Abstract 

The Marcus-Lushnikov process is a finite stochastic particle system in which each particle 
is entirely characterized by its mass. Each pair of particles with masses x and y merges into 
a single particle at a given rate K(x,y). We consider a strongly gelling kernel behaving as 
K(x, y) = x a y + xy a for some a G (0, 1]. In such a case, it is well-known that gelation occurs, 
that is, giant particles emerge. Then two possible models for hydrodynamic limits of the 
Marcus-Lushnikov process arise: the Smoluchowski equation, in which the giant particles are 
inert, and the Flory equation, in which the giant particles interact with finite ones. 
We show that, when using a suitable cut-off coagulation kernel in the Marcus-Lushnikov 
process and letting the number of particles increase to infinity, the possible limits solve either 
the Smoluchowski equation or the Flory equation. 

We also study the asymptotic behaviour of the largest particle in the Marcus-Lushnikov process 
without cut-off and show that there is only one giant particle. This single giant particle 
represents, asymptotically, the lost mass of the solution to the Flory equation. 

Keywords : Marcus-Lushnikov process, Smoluchowski's coagulation equation, Flory's model, gela- 
tion. 

MSC 2000 : 45K05, 60H30. 

1 Introduction 

We investigate the connection between a stochastic coalescence model, the Marcus-Lushnikov pro- 
cess, and two deterministic coagulation equations, the Smoluchowski and Flory equations. Recall 
that the Marcus-Lushnikov process [71 [5] is a finite stochastic system of coalescing particles while 
the Smoluchowski and Flory equations describe the evolution of the concentration c(t, x) of par- 
ticles of mass x £ (0,oo) at time t > in an infinite system of coalescing particles. Both models 
depend on a coagulation kernel K(x,y) describing the likeliness that two particles with respective 
masses x and y coalesce. When K increases sufficiently rapidly for large values of x and y, a sin- 
gular phenomenon known as gelation occurs: giant particles (that is, particles with infinite mass) 
appear in finite time (see Jeon [6], Escobedo-Mischler-Perthame [3]). There is however a strong 
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difference between the Smoluchowski and Flory equations: for the former, the giant particles are 
inert, while for the latter, the giant particles interact with the finite particles. 
When K(x,y)/y — > as y — > oo for all x G (0, oo), it is by now well-known that the Marcus- 
Lushnikov process converges to the solution of the Smoluchowski equation when the number of 
particles increases to infinity (see, e.g., Jeon jfj] and Norris [5]). On the other hand, it has been 
shown in that, if K(x,y)/y — ► l(x) G (0, oo) as y — > oo for all x G (0,oo), then the Marcus- 
Lushnikov process converges to the solution of the Flory equation. 

Our aim in this paper is to study more precisely how this transition from the Smoluchowski 
equation to the Flory equation arises in the Marcus-Lushnikov process. For a coagulation kernel 
K of the form K(x,y) ~ xy a + x a y for some a G (0, 1], we consider a Marcus-Lushnikov process 
starting with n particles, with total mass m n , where coalescence between particles larger than 
some threshold mass a n is not allowed. We show that, in the limit of large n, m n and a n , this 
Marcus-Lushnikov process converges, up to extraction of a subsequence, either to the solution of 
the Flory equation or that of the Smoluchowski equation, according to the behaviour of a n /m n for 
large values of n. 

We also study the behaviour of the largest particles in the Marcus-Lushnikov process without 
cut-off, and show that, in some sense, the total lost mass of the Flory equation is represented by 
one giant particle in the Marcus-Lushnikov process. Aldous pQ proved other results about giant 
particles for some similar (but more restrictive) kernels. We in fact obtain a much more precise 
result about the size of the largest particle after gelation, but we are not able to extend to our 
class of kernels his result about the largest particle before gelation. 



2 Main result 

Throughout the paper, a coagulation kernel is a function K : (0, oo) 2 i— > [0, oo) such that K(x, y) = 
K(y,x) for all (x,y) G (0, oo) 2 . We denote by M.f the set of non-negative finite measures on 
(0, oo). Let us first recall the definition of the Marcus-Lushnikov process. 

Definition 2.1 Consider a coagulation kernel K, and an initial state /xq = m 1 S Xf , with 

n > 1, (jci, x n ) G (0, oo)' 1 and m = x\ + ... + x n . A cadlag M~t -valued Markov process (nt)t>o 
is a Marcus-Lushnikov process associated with the pair (K,/j,q) if it a.s. takes its values in 

S(n,m):=l— ^S Vi , 1 < k < n, (y;)i<;< fe G (0, oo) fe , J^yi = ml (2.1) 

I m i=l i=l J 

and its generator is given by 

L K ^^{p) =]T{^ + m- 1 (5 yt+yi - 6 yt - <y ] - VM} (2-2) 

for all measurable functions ip : A4~^ i— » K and all states /i = mT 1 J3 i=1 S yi G S(n, m). 

This process is known to be well-defined and unique, without any assumption on K, see, e.g., 
Aldous O Section 4] or Norris Section 4] . 



2 



We now describe the Smoluchowski and Flory coagulation equations and first introduce the class of 
coagulation kernels to be considered in the sequel. As already mentioned, we will deal with kernels 
of the form K(x, y) ~ x a y + xy a for some a £ (0, 1]. More precisely, we assume the following: 

Assumption (A a ): The coagulation kernel K is continuous on (0, oo) 2 and there are a £ (0,1], 
I £ C((0, oo)), and positive real numbers < c < C < oo such that 

lim K(x,y)/y = l(x) , 

and 

c [x a y + xy a ) < K{x,y) < C {x a y + xy a ) and cx a <l(x) <C x a (2.3) 
for all (x, y) £ (0, oo) 2 . 

For such coagulation kernels, weak solutions to the Smoluchowski and Flory coagulation equations 
are then defined as follows: 

Definition 2.2 Consider a coagulation kernel K satisfying (A a ) for some a £ (0, 1] and /j,q £ 
such that (fio(dx), 1 + x) < oo. For (f> : (0, oo) i— » R, set 

A(j>{x,y):=(j } (x + y)-^{x)-^(y). (2.4) 

A family (fJ,t)t>0 C such that t i— * (fit{dx),x) and t ^ (/it(dx),l) are non-increasing is a 

solution to: 

(i) the Smoluchowski equation (S) if 

1 f 1 

(Mt,W = (^o,0) + 2y d J -s{dx)fi s (dy) 7 K(x,y)A(f>(x,y)) ds (2.5) 

for all <p £ C c ([0, oo)) and t > 0; 

(ii) the Flory equation (F) if for 

1 f* 

(lH,<P) = <A*o,0) + 2/ (l J 's(dx)fi s {dy) 7 K(x,y)A(f>(x,y)) ds 
t 

(fj, s (dx), <f>(x)l(x)) (l4)(dx) - fJ, s (dx), x) ds (2.6) 

for all 4> S C c ([0,oo)) and t > 0. Here and below, C c ([0,oo)) denotes the space of continuous 
functions with compact support in [0, oo). 

Note that the assumptions on K, (/it)t>o and <j> ensure that all the terms in (|2.5p and (|2.6[) make 
sense. 

Applying (|2.5[) (or (|2.6p ) with ^i(x) = x (which does not belong to C c ([0,oo))) would clearly give 
A</> = 0. Hence the total mass (/X{(rfx),x) is a priori constant as time evolves. However, for 
coagulation kernels satisfying (A a ), the gelation phenomenon (that is, the loss of mass in finite 
time, or, equivalently, the appearance of particles with infinite mass) is known to occur [31 [5], 
which we recall now, together with other properties. 
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Proposition 2.3 Consider a coagulation kernel K satisfying (A a ) for some a £ (0, 1] and /io £ 
such that (fio(dx), 1) < oo and (no(dx), x) = 1. For any solution (fit)t>Q of the Smoluchowski 
or Flory equation, the gelation time 

T gel := inf{f > : (n t (dx),x) < (fi (dx),x)} (2.7) 

is finite with the following upper estimate (here c is defined in (A a ) ) 

(Mdx)^ 1 -*) 
gel ~ (l-2-»)c ' 

If (fit)t>o solves the Flory equation, then t i— > (fi t (dx),x) is continuous and strictly decreasing on 

(T geh OO), 

/ 

for all e > 0. 



lim (n t (dx),x) =0 and I (/j, s (dx), x 1+Q ) ds < oo 



The proof that gelation occurs is easier under (A a ) than the general proof of Escobedo-Mischler- 
Perthame [3], and we will sketch it in the next section. This result expresses that particles with 
infinite mass appear in finite time. Observe next that equations (S) and (F) do not differ until 
gelation. The additional term in equation (F ) represents the loss of finite particles with mass x, 
proportionally to l(x) and to the mass of the giant particles (jj,o(dx) — n s (dx),x). 
Note that we are not able, and this is a well-known open problem, to show that t \— > (fj, t (dx),x) is 
continuous at t = T ge i. 

We finally consider a converging sequence of initial data. 

Assumption (J): For each n £ N \ {0}, we are given fift — m" 1 Y^i=i ^>x n for some (x™, x") £ 
(0, oo)™ and m n = + We assume that there exists /io € such that (/j,o(dx),x) = 1 and 

lim„ (/j,q (dx), 4>) — (no(dx),(f>) for all <fi £ Cf,([0, oo)), C&([0, oo)) denoting the space of continuous 
and bounded functions on [0, oo). In addition, 

lim sup (^(tfc),!^]) =0. 

We will actually not use explicitly all the assumptions in (J) and (A a ): some are just needed to 
apply the results of [5]. We now state a compactness result which follows from [5]. 

Proposition 2.4 Consider a coagulation kernel K satisfying (A a ) for some a £ (0,1] and a se- 
quence of initial conditions (hq)u>i satisfying (I). For each a > and n > 1, we put K a := 
K 1 

(o,o]x(o,a] an d denote by (fi^ a ) t >o the Marcus- Lushnikov process associated with the pair 
(Ka,^). The family {(^™' a )t>o}n>i,a>o is tight in B([0, oo), .M ^ ), endowed with the Skorokhod 
topology associated with the vague topology on M.^ . 

This proposition is proved in [5l Theorem 2.3-i] (with the choice of the subadditive function <p(x) = 
V2C(1 + x), for which K(x,y) < 4>(x)4>(y)). Actually, it is stated in [5] without the dependence 
on a, but the extension is straightforward. 

Notice here that, if a > m n , the Marcus-Lushnikov process (//™' a )t>o reduces to the standard 
Marcus-Lushikov process associated with (K,[Iq ). 

We may finally state our main results. Recall that we assume the total mass of the system to be 
initially (/xq (dx) , x) = 1. 
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Theorem 2.5 Consider a coagulation kernel K satisfying (A a ) for some a € (0, 1] and a sequence 
of initial conditions (/Uo)n>i satisfying (J). Consider also a sequence (a n ) n >i of positive real 
numbers such that lim„ a n = oo. For each n > 1, Zei (/i"' a ")t>o fre i/ie Marcus- Lushnikov process 
associated with the pair (K an , where K an := -K" l(o,a n ]x(o,a„] and consider the weak limit (/it) t>o 
in D([0, oo), y) o/ a subsequence |(/i™'" a ™' c )t>o| ■ Then (fi t ) t >a belongs a.s. to C([0, oo), Al y) 
and enjoys the following properties: 

1. Assume that a n = in n . 

(i) Then (nt)t>o solves a.s. the Flory equation with coagulation kernel K and initial condition 
Mo- 

(ii) Denote by M"(t) > M%(£) > ... £/ie ordered sizes of the particles in the Marcus-Lushnikov 
process (M™' m ™)t>o> a* 7 ^ define the (a priori random) gelation timeT ge i of(/j,t)t>o as in \2. 7| j. 
TTien /or all rj > and (3 > 0, 



lim £ 



sup 

te[T 9ei +r ): oo) 



M™ k (t) 



- (1 - (fi t (dx),x)) 



= 0, 



lim lim sup P 

6— »oo fc^oo 



f° — ^J^( S )l M (Mf( S ))(l S >/? 



(2.8) 
(2.9) 



Furthermore, there is a positive constant L depending only on K such that, for all rj > and 
b>l, 



lim E 



M? k (t) 
sup — 

te[o,r 9 ei-i7] TO ™fc 



0, 



lim sup E 

k — >oc 



jJ-s 



'{dx),x\[ b}00 ){x)) ds 



< 



(2.10) 
(2.11) 



2. If a n /m n —> as n —> oo, then (fit)t>o solves a.s. the Smoluchowski equation with coagulation 
kernel K and initial condition fio . 

3. If a n /m n — > 7 S (0,1) as n — > oo, then {(J't)te[o,T 1 ) solves a.s. the Flory equation with 
coagulation kernel K and initial condition /j,q where 

Ti := inf{t > : 1 - (fit(dx), x) > 7}. (2.12) 

Point l-(i) is proved in [5J Theorem 2.3-ii]. Remark that (|2.10[) is almost obvious while (|2.1ip gives 
an estimate on the tail of the mass distribution before gelation. The most interesting estimate 
is of course (|2.8I) which shows that, for t > T ge i, the largest particle in the Marcus-Lushnikov 
process without cut-off occupates a positive fraction of the total mass of the system with a precise 
asymptotic. Finally, (|2.9p shows that, in some sense, there is only one giant particle after gelation: 
the other particles are rather small. Other results about the largest particles for the kernel 

K< \ = 2(^) 1+ ° 
{XlV) {x + y^+a-x^-y^ 1 
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which satisfies (A a ), were obtained by Aldous [T]. He however did not show that, after gelation, 
the size of the largest particle is of order em n . 

Point 2 seems to be new, and quite interesting. Indeed, we allow arbitrary cut-off sequences (a n ) 
which increase more slowly than (m n ). 

Finally, Point 3 can be explained in the following way: assume that a n = ~/m n for all n > 1 and 
some 7 € (0,1) and that there is only one giant particle in (/i"' m ™) t > . In that situation, we 
then clearly have (M™' 7m ")te[o,T 1 "] = (/ i t l ' m ")te[o,T™] I where T™ is the first time at which the giant 
particle has a size greater than 7rn n , i.e., it occupates a fraction 7 of the total mass of the system. 
Thus, (£i"' 7m ™)te[o,T 1 "] should converge to {^t)te[o,T!], where \i solves the Flory equation, and T\ 
is the first time for which the giant particle occupates a fraction 7 of the total mass in the Flory 
model. 

The proof of Theorem 12.51 is given in Section [4] after establishing some properties of solutions to 
the Smoluchowski and Flory coagulation equations in the next section. The final section of the 
paper is devoted to numerical illustrations. 



3 Properties of solutions to (S) and (F) 

Throughout this section, K is a coagulation kernel satisfying (A a ) for some a £ (0,1] and Hq 
belongs to M.^ with total mass (no(dx) , x) = 1. 

Proof of Proposition 12.31 Let (/^t)t>o be a solution to the Smoluchowski equation (S) or the 
Flory equation (F), and define T ge i <E (0, 00] by (|2.7p . Classical approximation arguments allow 
us to use (|2.5p and (|2.6[) with <j>(x) = x 1 ^". Indeed, it suffices to approximate by a sequence of 
functions in C c ([0, 00)) and to pass to the limit, using the first inequality in 

mints:, y)^ a > x^ a + y 1 - - (x + yf- a > (2 - 2 1 "") mm(a;, y) 1 ' 01 , (3.1) 

which warrants that K(x,y)\A(j)(x,y)\ < 2Cxy by (A a ). We deduce from (|2.5p . (|2.6p . and the 
second inequality in (|3.ip that, for all t > 0, 

(^{dx),x 1 - a ) < (fioidx)^ 1 - 01 ) 

2 J (^s(dx)ii s (dy),K(x 1 y)mm(x,y) 1 ~ a ) ds. 

By virtue of (|2.3p . K(x,y) min(x, y) 1 "" > cxy, whence 

(l-2- a )c [ (ii s (dx),x) 2 ds < (noidx)^ 1 - 01 ) (3.2) 

for all t > 0. Since (fj, s (dx), x) — (fj,o(dx), x) = 1 for all s € [0, T ge i), we realize that T ge i has to be 
finite for (|3.2p to hold true. A further consequence of (|3. 2|) is that 

(1 - 2- a )cT gel = (1 - 2~ a )c [ (v s {dx),x) 2 ds < {^{dx)^ 1 -*) , 

Jo 

whence 

9el ~ (l-2-«)c ' 
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It also follows from (|3.2p that t i — ► ((At(dx), x) belongs to L 2 (0, oo) which, together with the 
monotonicity and non-negativity of i 1 — > (fi t {dx), x) implies that (fi t (dx),x) — ► as t — ► oo. 

We now assume that (/X{)t>o solves the Flory equation (F) and prove that, for all e > 0, 

(^t(da;),a; 1+Q ) eft < oo. (3.3) 

T gci +e 

To do so, we apply (|2.6p with the choice 4>a{%) = niin(x, A) for some positive real number A. Since 
A<p a is non-positive, we get 

(Hu<Pa) < (ho,4>a) - / (ju s ((ia;),min(a;, A)l(x)) (1 - {fx s (dx),x)) ds. 
Jo 

Since 1 = (/j,a(dx), x) > (n s (dx),x) we may let A — ► oo and f — > oo in the above inequality and use 
the Fatou lemma to deduce that 



(fi s (dx),xl(x)) (1 - (fi 3 (dx),x)) ds<l. (3.4) 
Let e > 0. On the one hand, putting 

it follows from the definition (|2.7p of T se j that <5 E > 0. On the other hand, xl(x) > cx 1+a by (A a ). 
We therefore infer from (13.411 that 



(^(dx),x 1+Q )ds< 4~, 
whence p.3p . 

We now check that t i— ► (fi t (dx), x) is continuous on (T se ;,oo). Using once more (|2.6p with the 
choice ^(x) = min(x, A), we obtain for T ge i < s < t 



1 



(fit- fj, s ,(p A ) = -J dj, T (dx)ii T (dy),K(x,y)A(pA(x,y))dT 
t 

(^ T (dx),tp A (x)l(x)} (1- (fj, T (dx),x))dr. 

Clearly A<Pa(x, y) — > as A — ► oo for all (x, y) G (0,oo) 2 while (A a ) warrants that 
K(x, 2/)|A^(x, y)\ < C(x a y + xy a ) min(x, y) < C(x 1+Q y + xy 1+a ). 
Using (|3.3p and the Lebesgue dominated convergence theorem, we obtain 

(Ht{dx) — (j. s (dx), x) = — J (fj, T (dx), xl(x)) (1 — {(j. T (dx), x)) dr . (3-5) 

Using again (|3.3[) and that | (fj, T (dx), xl(x)) (1 — (fi T (dx), x)) \ < C (fj, T (dx), x 1+a ) by (A a ), we 
conclude that 1 1 — ► (fi t {dx), x) is continuous on (T ge i, oo). 
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It remains to check that t i — > (fit(dx),x) is strictly decreasing on t € (T ge i, oo). According to 
(|3.5p this is true as long as /t r 7^ for T ge i < s < r < t: it thus suffices to show that /it for 
all i > 0. For that purpose, we take <f>(x) = x1(o.a]( x ) hi (BUI) where A > is chosen so that 
(no(dx), x1(o j a] (#)) > (such an A always exists as (no(dx),x) — 1). Thanks to (|2.3p . we get 

> -C {^{dx^tidy), (x a+1 y + x 2 y a )l [0A] (x)) - C (» t (dx), x 1+a l {0>A] (x)) 

> -C(A a + A) (/it, <£) ( Mt (dx), a; + z Q ) - CM" ( Mt , <t>) ■ (3.6) 

Since t 1 — ► (fi t (dx),x) and £ 1 — ► (fj, t (dx),x a } are non-increasing and (/io(dx), x + x a ) < 00, we 
conclude that 

-Tj_ <t>) > ~Ca (1 + (fJ, (dx),x a )) (fi t , 4>) > -C A>IM > ((H,4>) 

for all t > for some constant Ca,h > 0. Consequently, (fit, </>) > for all t > as the choice of 
A warrants that (/to, </>) > 0, and the proof of Proposition 1 2 . 31 is complete. □ 

Next, as a preliminary step towards the proof of Theorem 12.51 Point l-(ii), we show that solutions 
to the Smoluchowski and Flory coagulation equations do not coincide after the gelation time. 

Corollary 3.1 Let (nt)t>o ond (vt)t>o be solutions to the Smoluchowski equation (S) and the 
Flory equation (F), respectively, (with the same coagulation kernel K and initial condition /j,q), 
and assume further their respective gelation times coincide, that is, 

T ge i := inf{f > : (fi t (dx),x) < (fx (dx),x)} = inf{t > : (v t (dx),x) < (fj, (dx), x) } . 

Then, for each e > 0, there exists s E S (T ge i,T ge i + e) such that fj, Se 7^ v Se . 

Proof. Consider e > 0. 

Either t 1 — ► (/j, t (dx), x 1+a ) does not belong to L 1 (T ge i + (e/2),T ge i + e) and /it cannot coincide 
with v t on (T gel + (e/2), T gei + e) since t 1 — ► (v t {dx), x 1+a ) belongs to L 1 (T gel + (e/2), T gel + e) 
by Proposition ^. 31 

Or t 1 — > (nt(dx), x 1+a ) belongs to L 1 (T ge i + (e/2),T ge i + e) and it is not difficult to check that 
this property and (|2.5p entail that (/j, t (dx), x) = (^T gel +(e/2)(dx), x) for t £ [T ge i + (e/2), T ge i + e]: 
indeed, take <Pa(x) = min(x, A) in (12. 5|) and pass to the limit as A — > 00 using that A<p A (x, y) — > 
and the time integrability of t 1 — > ( x /it(dx),x 1+a '). Owing to the strict monotonicity of t 1 — > 
(v t (dx),x) established in Proposition ^. 31 the previous property of /it excludes that /t t = f t for all 
t 6 (r ge ; + (e/2), T ge i + e) and completes the proof of Corollary 13. II □ 



4 Proof of the main results 

We fix a coagulation kernel K satisfying (A a ) for some a S (0, 1] and a sequence of initial data 
(Mo)»>i satisfying (7). Next, for a > and n > 1, we put K a (x,y) := if (#, y)l(o,o]x(o,o] an d 
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denote by (/x"' a ) t >o the Marcus-Lushnikov process associated with the pair (i^/ip). According 
to Definition 12.11 we may write 



— £<Wr°M with Mr(t)>M^ a (t)>MS' a (t)>... (4.1) 



for all t > 0, n > 1, and a > 0. 

We next recall that the space O([0, oo), ) is endowed with the Skorokhod topology associated 
with the vague convergence topology on (see Ethier-Kurtz [lj for further information), and 
denote by d a distance on M. ^ metrizing the vague convergence topology. 

Marcus-Lushnikov processes have some martingale properties, which are immediately obtained 
from (|2"T2")1 , see also Section 4]. 

Lemma 4.1 For all cf> £ L^ c (0, oo) and t > 0, we have 

( M r^) = ^ i »+or(0) 

1 r Y, K a{M^ a {s),M^ a {s))^{M^ a {s),M^ a { S )) ds 



2 ™n Jo 



(l$,<i>)+0? a (</>) + - / (jJp a (dx)tf a {dy),K a (x,y)A<l>(x,y)) ds 



1 



(f^< a (dx),K a (x ) x)A(f)(x,x)) ds 



(4.2) 



2m n jo 

where Aip is defined in \2.4\ l, and O n,a ((j)) is a martingale starting from with (predictable) 
quadratic variation 



(O n ' a (<t>)) t 



1 



2j7l n Jo 

1 



H? a (dx)(iZ> a (dy),K a {x,v) [A^{x,y)Y) ds 



^ a (dx),K a (x,x)[Acf > {x,x)} 2 ) ds 



Jo 

Furthermore, if <p : (0, oo) — > R is a subadditive function, that is, (f>(x + y) < 4>{x) + <f>(y) for 
(x,y) € (OjOo) 2 , then 1 1— > (/i™' a ,0) is a.s. a non-increasing function. 

We carry on with some easy facts. 

Lemma 4.2 Lei (a n ) n >i oe a sequence of positive numbers. Then any weak limit (^t)t>o of 
the sequence {(Ai"' a ")t>o}„> 1 belongs a.s. to C([Q, oo), M~f ), and both t \— > (fit(dx),x) and t i— > 
(fit(dx),l) are a.s. non-increasing functions. Furthermore, 

(4.3) 



sup sup (fi"' a " (dx), 1 + x) = k := sup (fifi (dx), 1 + x) < oo, 

n>l t>0 n 



and for all <f> £ C c ([0, oo)) and T > 0, 



lim E 

n — >oo 



lim E 



sup 

*e[o,T] 



2m r , 



(jj^' an (dx), K an (x, x)A<j>(x, x)) ds 



sup (or-wr 

te[o,r] 



0. 



(4.4) 
(4.5) 



9 



Proof. First, if <j> £ C(,([0, oo)), the jumps of (/x™' a ",(/>) are of the form m" 1 A</(>(x, y) and clearly 
converge to zero as n — » oo since m n — ► oo. Hence any weak limit (nt)t>o belongs to C([0, oo), .Ml") 
a.s. 

Consider next a family (Afc)b>o of continuous non-increasing functions on (0, oo) such that Xb{x) — 
1 for x < b and Xt,(x) = for x > 6+1, and a non-negative subadditive function 0. Then, on 
the one hand, tftXt, is also subadditive and Lemma |4. 1 1 ensures that t 1 — ► (m"' ") '^fr) is a - s - non- 
increasing for all n > 1. On the other hand, since 4>Xb € C c ([0, oo)), it follows from the definition 

of (fJrt)t>o that there is a subsequence (rik)k>i, nk — * °°, such that < ( ( /x? fc ' a " fc , <f>Xb ) ) > 

converges in law towards ((/it, </><-tb)) t >o for each fixed b > as k — * oo. Therefore, 1 1 — ► (/it, 
is a.s. non-increasing for each b > 0. Since ((/Xt, 0A'b)) f)>o converges to (fHifi) as 6 — > oo for each 
t > 0, we conclude that £ i— » (/it,0) is a.s. non-increasing. Applying this result to </>(x) = 1 and 
</>(x) = we obtain that both t 1— > (fi t (dx), x) and 1 1— > (fi t (dx), 1) are a.s. non-increasing functions 
of time. 

Next, since a; i— > 1 + a; is subadditive, Lemma 14.11 implies that we have a.s. (/i"' a " (dx), 1 + x) < 
(fj,Q(dx), 1 + x) for n > 1 and t > 0, and (/ip (dx), 1 + a;) is bounded uniformly with respect to n 
by assumption (7). 

Consider finally <j> G C c ([0, oo)) with support included in [0, R] for some R > 0. By (|2.3[) . 

\K an (x,x)A(t>(x,x)\ < 6C||0||ioo whence 



| (^' a " (di), i^a„ (as, x)Act>(x, x)) | < 6Cre|| 



| L oo i? 1+Q a.s. 



by (|4.3p . from which (|4.4p readily follows since m n — ► oo. By a similar argument, we establish that 
E [(O n,CLn ((/))) t ] — ► as rn oo, which implies (|4.5[) by Doob's inequality. □ 



We now prove a fundamental estimate which provides a control on the large masses contained in 

n,a 

/V • 

Lemma 4.3 There exists a positive real number L depending only on c and a in (A a ) such that 

L 



E 



i^Mf< Q ( S )M;' Q ( S )l [M (Mf'») (M;<») ds 



< 



(4.6) 



for all n > 1, a > 0, and b € (0, a), the M"' a kmj defined in {^.1 



Proof. To prove this estimate, we use (|4.2[) with <j)(x) = x 1_Q min (x, b) a for some 6 € (0, a). We 
first notice that (/j,q,4>) < (Ho (dx) , x) = 1 and (//™' a ,0) > for all t > and n > 1. In addition, 
is subadditive so that A<j)(x, y) is always non-positive and we infer from (|2.3p and (|3.ip that 

/f a (a;,y)A0( a; ,y) < -(2-2 1 - Q )c6 Q xyl [M (x)l [M (y) 

for (x, y) € (0,oo) 2 . Taking expectations in (14. 2p and using the above inequalities, we obtain 



< 1 - — E 

Li 



-L^Mr( S )Mf a ( S )i [M] (Mr( S ))i M (Mf a ( S )) d S 
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for all t > 0, with 1/L := c(l — 2 Q ). We conclude the proof by letting i — > oo in the previous 
inequality. □ 

We now turn to the proof of Theorem 12.51 and first recall that Point l-(i) is included in [5J Theo- 
rem 2.3-ii] as (/x™' m ")t>o is the standard Marcus-Lushnikov process associated with (K,[1q). 

Proof of Point 2 of Theorem 12.51 Let (a ra ) n >i be a sequence of positive real numbers sat- 
isfying a n — > oo and a n /m n — ► as n — ► oo. We consider the limit (nt)t>o of a subsequence 
i (^J**" "* ) t > I m f-jjg sense that a.s. 

lim sup d( f i" k ' a " k ,tH)=Q for all T > , (4.7) 

the existence of such a limit being guaranteed by Proposition 12.41 Lemma 14.21 (which ensures the 
time continuity of fit) and the Skorokhod representation theorem. We now aim at showing that 
{Ht)t>o solves a.s. the Smoluchowski equation (5) and proceed in two steps. 
Step 1. We first deduce from Lemma FOl that 



E 



T 2 

(n™' an {dx),x\[ b ^ an ](x)) dt 



<( 7^+^—1 (4-8) 



b a in 

for all b > 0, n > 1, and T > 0. Indeed, we have a.s., for all t > 0, 

= (/ i r' Q "( da; )Mr Q "(^) J a; y a [ft,a„](a;)l[b,a„](y)) - — <M™' a "((ix),X 2 l [i) , aii j(x)) 

> (/i™' Q "(da;),x]l [f , ia7i ](a;)) 2 - — {^ an {dx),x) (4.9) 

> (^™' Q "(dx),a;lL [b ri] (a;)) 2 - — , 

hence (I4.8[) after integrating over (0,T), taking expectation, and using Lemma 14.31 (with a — a n ). 

Step 2. By Lemma I4.2[ we already know that t \— > (fit(dx),x) and t i— > (fi t (dx), 1) are a.s. non- 
increasing functions. Consider now </> G C c ([0,oo)). The convergence (|4.7|) and the assumption 

(I) ensure that ^/i™'" a ™ fc j ^ — > (/■«*> <^>) a.s. for all t > and (/iQ fc ,(/>) — ► (fio,(j)} as fc — > oo. 
Recalling (|4.2[) . (|4.4[) . and (|4.5[) . we realize that (fit)t>o solves (|2.5[) provided that we check that 
-Bfc(t) — * B(t) (for instance in i 1 ) as k — ► oo for alH > 0, where 

B fc (t) := ^ ^T' a " k (dx)fj 7 s k ' a ' 1 ' k {dy) 7 K aiik {x,y)A<j)(x,y)^ ds, 

B(t) := / (fi s (dx)fi s (dy),K(x,y)A(f>(x,y)) ds. 
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For that purpose, we consider a family (Xb)b>a of continuous non-increasing functions on [0, oo) 
such that X b (x) — 1 for x € (0, 6] and X b {x) — for x G [6+1, oo), and put 

B k (t,b) := J ( f i: k ' a ^(dx) f iT' a " k (dy),K ank (x,y)Acj > (x,y)Xb(x)Xb(y)) ds, 

B(t,b) := f (n s (dx)n s (dy),K(x,y)A^(x,y)X b {x)X b (y)) ds. 



On the one hand, it follows from (A a ), the boundedness of 0, the bounds x a < 1 + x and 

sup (nt(dx), 1 + x) = (no(dx), 1 + x) , 

t>0 

and the Lebesgue dominated convergence theorem that 

lim E[\B(t,b) - B(t)\] = 0. 

b — ^oo 



(4.10) 



On the other hand, for each 6 G (0, oo), we have K a (x, y)Xb{x)Xb{y) — K(x, y)Xb(x)Xb(y) for all 
(x,y) £ (0,oo) 2 as soon as 6 + 1 < a„ fc , the latter being true for k sufficiently large. Consequently, 
since (x, y) i — ► K(x,y)Xb(x)Xb(y)A(f>(x,y) belongs to C c ([0,oo) 2 ), the convergence f|4.7[) entails 
that Bk(t, 6) — ► B(t, 6) for alH > a.s. as k — > oo. Thanks to (A a ) and (|4. 3|) we may apply the 
Lebesgue dominated convergence theorem to obtain that 



lim E[\B k (t, 6) - B(t, 6)|] = for each 6 > . 



k— »oo 



(4.11) 



Finally, owing to (|2.3p . we have for k sufficiently large (such that a nk > 6) 



E[\B k {t,b)-B k (t)\ 

rt 

<SM\\ l »e' 



nh,a„. 
Us 



< 6C\\cf>\\ L °oE 

< 6C\\<P\\ L ~E 
+6C\\<t>\\ L ~E 

< QC\\<j)\\ L -E 
^U\\l- e 



(dx)nT' a ^ (dy),K ank (x,y) (1 - X b {x)X b {y))) ds 
(ji n s k ' ank {dx)fxT' ank (dy),x a yl i0 , ank] (x)l {0yank] (y) (1 - X b (x)X b (y))^ ds 
^ n s k ' an " (dx)fiT' a71k (dy),x a yl {0 ^ nk] (x)l [b ,a nk ](y)^ ds 
fj^ k ' a " k (dx)(i™ k ' ank (dy), a; Q ?/l[^ Q „j(x)l(o,Q„j(y)^ ds 
^s k ' a,lk (dx),(l +a;)l( , a „ fe ](x)\ ^" fc ' a " fe (dy),yl[ biank ](y)\ ds 
[iT' a " k (dx),xl [btai ](x)) (^r' Q " fc (d?/),y]l(o,a„,](y)) ds 
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We then infer from (|4.3[) . the Cauchy-Schwarz inequality, and (|4. 8[) that, for b > 1, 
S[|B fc (t,6)-B fc (t)|] < 12/«C||0|| i oo J E ^ ^""'^(rfa;)^!^]^)) ds 

< UkCUWlS^E J ^T ,ank {dx),xTL [btank] {x)y ds 

\0 m n k / 



1/2 



1/2 



Since a n /m n — + as n — > oo, we may first let A; — * oo and then 6 — » oo in the previous inequality 
to conclude that 

lim limsup.E[|B fc (t,&) - B k (t)\] = 0. (4.12) 



Combining (|4T0|) . (|4~TT]) . and (f47T2|) ends the proof. 



We next complete the proof of Point 1 of Theorem 



□ 



Proof of Point of Theorem 12.51 Recall that we are in the situation where a n — m n , 

so that (^("' a ") t >o = (/i"' m ")t>o is the classical Marcus-Lushnikov process associated with (K,/Iq) 

for each n > 1. Let ([it)t>o be the limit of a subsequence <. (/i t : " k )t>o \ in the sense that a.s. 

I J &>i 



lira sup d(^ k ' m " k , /^) = for all T > . 



(4.13) 



the existence of such a limit following by the same arguments as (|4.7j) . 

We already know from [5] that (nt)t>o solves a.s. the Flory equation. We define the (a priori 
random) gelling time T ge i of (nt)t>o by (|2.7p and write 



Mi 



" = with M "W > M 2 "(t) > M 3 *(t) > 



(4.14) 



for all t > and n > 1. 

As before we denote by (Af,)b>o a family of continuous non-increasing functions such that Xb(x) = 1 
for x £ [0,6] and X b (x) = for x > b + 1. We start with the proof of (|2.10[) which is almost 
immediate. Since 1 i— > M"(t)/m n is a.s. non-decreasing and bounded by 1, it suffices to check that 
P{M™ k (T ge i - rf) > 5m„ k ] — ► as k — ► oo for all 77 > and <5 > 0. For that purpose, fix b > 0. 
Since 1 - Af 6 > 0, it follows from P~T4"|) that 



(daj),a:(l- X b (x))} > 
(dx),xXb(x)) > 



n k ,m n 



M 



M 



nk,m„ k 



(*). 



L [b+l,i 



3) (t)) 
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For k large enough we have Sm nk > b + 1 and thus 



E 



1 ( ^Tgli-v ^ x )i xX b(x) 



> E 



M^{T gel - n ) 



[8m n ,00) 



(m,— (T gel -r,)) 



> 6P[M?*(T gel -ri)>6T 



Now, thanks to the compactness of the support of X b and (|4.13|) . the sequence 



k>l 



converges a.s. to (fXT ge i-r](dx),xXi ) (x)) and is bounded by 
above inequality to obtain 



We may then let k — > 00 in the 



E [1 - (ti Tgsl - v (dx),xX b (x))] >SlimsupP[M? k (T gel -ri)>5m nk } 



k — ►c 



Next, owing to the definition of T ge i, the sequence ^fix ge i-ri(dx), xXb(x))) b converges towards 
(/j,T g ei-r)(dx), x) = 1 as b — > 00 and is bounded by 1. Passing to the limit as b — > 00 in the previous 
inequality entails that P [M™ k (T ge i — rf) > Sm nh ] — > as k — > 00, which is the claimed result. 
The limit (f^TTU]) then follows. 



00, we have a Hk > b 



We now turn to the proof of (|2.11|) . Let b > 0. Since a n = m n — > 00 as n 
for k large enough and it follows from Lemma 14.31 (with a = a n ) by an argument similar to (|4.9[) 
that (recall that all the particles represented in fi t ' " fc are smaller than m nk by construction, see 
Definition O) 



E 



<—+E 



n k ,m n 



(dx),xl [by00) (x)\ dt 



-gel J. 



< 



< 



L 
L 



E 



E 



Tgel ]\/f n k 



M" fc (t) 



n k ,m n 

Mt 



(dx),x 2 H [by0o) (x)J dt 
(dx),x) dt 



M" k (t) 



dt 



(4.15) 



Since M{ lk (t) < m„ fc and T ge i is a bounded random variable by Proposition ^. 31 we easily deduce 
from (|2~TUj) that E j^ gel (M? k (t)/m„ k ) dil — ► as fc -> 00. Thus [gXrj follows from (|435]) . 

We next establish (|2.8[) and (|2.9[) and split the proof into five steps. In the first two steps we show 
that, for t > T ge i, at least one particle has a size of order Sm n for some 5 > 0. Since such a particle 
is very attractive, we deduce in Step 3 that no other large particle can exist and obtain (|2.9p . We 
then conclude in the last two steps that, for t > T ge i, this single giant particle is solely responsible 
for the loss of mass and obtain (12.81). 
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Step 1. Let (a„)„>i be any sequence of positive numbers such that a n /m n — > as n — > oo. The 
aim of this step is to show that 

lim P[M?* {T gel + e) > a nk ] = 1 for all s > . (4.16) 

k—>OG 

For that purpose, we introduce the stopping time 

n :=M{t >0 : M**(t) > a„ J , 

and notice that we may assume that a n — > oo as n — > oo without loss of generality. Owing to the 
time monotonicity of M™, it suffices to prove that P[rk > T ffe z + e] — > as fc — > oo for all e > to 
establish (|4.16[) . Assume thus for contradiction that there is e > such that S :— limsup fe P[rk > 
Tgei + s] > 0. Then, on the one hand, we have P[ti] > S > with £1 := {limsup fe Tk > T ge i + e}. 
On the other hand, for each fc > 1, it is clearly possible to couple the Marcus-Lushnikov processes 
(/x t ' " fc )t>o and (p, t ' " fc )t>o in such a way that they coincide on [0, rfc). Hence, a.s. on Cl, up to 
extraction of a subsequence, 



lim sup d(/j, t 

k ^°° [0,T se! +e/2] 



fit) = lim sup d(fi' t 

k ^°° [0,T g<!l +e/2 



ih) = o. 



By Theorem 12.51 Points l-(i) and 2, we deduce that the limit (fJ>t)t>0 solves simultaneously the 
Flory and Smoluchowski equations on [0,T se i + e/2) with positive probability, which contradicts 
Corollary O 

Step 2. We now deduce from Step 1 that 



lim lim inf P {M? h {T get + e) > Sm rik ] = 1 for all e > . 

S—>Q k 



(4.17) 



Assume for contradiction that there is e > for which (|4.17p fails to be true. Then there exists 
7 G [0, 1) such that lim inf k P [M™ k (T ge t + e) > 6m nk ] < 7 for all 5 > 0. We may thus find a strictly 

increasing sequence (fc/)/>i such that P M 1 k ' (T ge i + e) > m nk /I < 7 for every I > 1. We then 

1 /2 

put a nk — m nk Jl for I > 1 (and e.g. a n — m n if n £ {n^ : Z > 1}). Then a n /m n — > as 
n — >■ 00 and the assertion (|4.16p established in Step 1 warrants that P [M™ k (T ge i + e) > a nk ] — ► 1 



as fc — > 00. But P 



Mi fci (T gei + e) > a n 



< 7 < 1 for all I > 1, hence a contradiction. 



.Step 5. We are now in a position to prove (|2.9|) which somehow means that the other particles are 
small in the sense that 



lim lim sup P 




; (s,6) ds >/3 



= 



(4.18) 



for all P > with the notation 

X n (s,b) 



1 



First note that a.s., for all s > and n > 1, we have 
-L ]T M? (s)M?(s)l [b , mn] (M?(s)) a [b)r 



M ;( S )) > -AA£2 1[i)i0o) (Mf( s )) x»( s ,6). 
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since M"(s) < m n for all s > and i > 1. By Lemma |4"31 (with a = m n ), we obtain 



E 



Mf(s) 



l[6,oc) (Mfta)) X"(s,6) 



< 



(4.19) 



We next fix (3 > 0,r]>0, and b > 0. By (j4~T7|) there is 6 > such that liminf fe P[M™ k (T gel +(3/2) > 
<5w nfc ] > 1 — 77. Recalling that 1 i— > M™(i) is a.s. non-decreasing and A™(s, 6) < 1 for all s > a.s., 
we have for k sufficiently large such that Sm rik > b, 



X nk (s,b) ds)>(3 



< P 



T 9 ei+/3/2 



X nk (s,b) ds > (3/2 



<P[M^(T gel +/3/2)<Sm nk ] 

M^{T gel +P/2) 
'T gcl +fj/2 Sm nk 
<1-P [M" k {T gel + (3/2) > 5m nk ] 



-P 



1 [6;00) (M? k (T gel +/3/2))X n "(s,b) ds > (3/2 



-P 



M™ k (s) 

T 3c! +/3/2 m n k 



l [b ,oo) {M? k {s))X nk {s,b) ds) >(35/2 



2L 



<1-P [M" k (T gel + (3/2) > 5m nk ] + — , 

the last inequality being a consequence of (|4.19p . Letting k — > oo in the above inequality, we 
obtain, thanks to the choice of S, 



lim sup P 

k — >oc 



X nk (s,b) ds > (3 



<T) + 



2L 

b a (38 



Now, we first pass to the limit as b — > oo and then as r\ — > in the above inequality to obtain 
i.e. (EU). 

Step 4. Set 7 4 := 1 - (n Tgel +t(dx),x) and B fc (t) := M" fc (T geZ + ()/m„ t for i > and k > 1. Our 
aim in this step is to prove that 



lim E 

k — >oo 



/ \B k (t) - lt \ dt 
Jo 



for all T > . 



(4.20) 



As before let (Xb)b>o be a family of continuous non-increasing functions such that Xb{x) = 1 for 
x G [0, b] and X b {x) = for a; > b + 1. We then put 

7 t b := 1 - (Mr ee ,+t(da;),a:^6(a;)) and 7 t 6,fe := 1 - ^/^I+t* (dx), avft,^)^ 

for 6 > 0, > 1, and t > 0. On the one hand, we have a.s. that 7 t b — > 74 as 6 — > 00 for all £ > 0. 
Since |7t| + I7* | < 2, we deduce from the Lebesgue dominated convergence theorem that 



lim E 



l7t-7tl dt 



0. 



(4.21) 
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On the other hand, owing to the compactness of the support of X b , we infer from (|4.13[) that 
7 t ' — ► 7j a.s. for all b > and t > 0. As |7t | + |7 t ' '| < 2, we use again the Lebesgue dominated 
convergence theorem to obtain that 



lim E 



h b t-7t' k \dt 



for each b > . 



(4.22) 



But 7 t b,fc = A btk (t) + B k (t) - C b>k (t) a.s., where 
A 6)fe (t) := 



C b , k (t) 



J_ £ ^.m.,, + t) ^ _ ^ (M ;— , {Tgei 

<!^k±i) A . t(M r-- m .. +< ))<^. 



Clearly, 



lim E 

k — >oo 



C btk (t) dt 



= for each b > . 



while, since < A b<k (t) < X nk (T gel + 1, b) < 1 a.s. 



lim lim sup E 



A b , k (t) dt 



= 



(4.23) 



(4.24) 



by (|4.18p and the Lebesgue dominated convergence theorem. 

Now, since B k {t) - 7t = C b , k [t) - A bik (t) + ( 7t b ' fc - 7 b ) + ( 7t fc - 7t ) for b > 0, k > 1, and * > 0, it 
follows from (|4T22]) and (|4~2^1l that 



lim sup E 

k — *oo 













7*1 dt 


< lim sup E 


/ A b>k (t) dt 


+ E 


r \it-it\dt 




k — >oo 


Jo 




JO j 



for all b > 0. Letting b — > oo and using (f¥^T]l and fl~2"4"j) give gHQ]) . 

Siep 5. To complete the proof of (|2.8|) . it remains to show that, for all e > and 77 > 0, we have 



lim P 

k — *oo 



sup \B k (t) - 7(| > e 

te[?7,oo) 



(4.25) 



Indeed, (|4.25|) clearly implies (12. 8p by the Lebesgue dominated convergence theorem since B k (t) < 1 
and 7 t < 1 for fj > a.s. 

We thus fix £ > and 77 > 0. Since (nt)t>o solves the Flory equation (F) by Theorem 1 2 . 51 Point 1- 
(i), it follows from Proposition 12.31 that t 1 — > 74 is a.s. increasing and continuous on [77,00) and 
7t — > 1 as t — > 00 a.s. It is also straightforward to check that 1 1— > -Bfc(i) is a.s. non-decreasing on 
[77, 00) with -Bfc(t) — > 1 as rj — > 00 a.s. As a consequence of the a.s. monotonicity and boundedness 
of 1 1 — ► 74 and i?^ (t) , we have for t >T 

\B k {t)- lt \ = max{B fc (i)-7T + 7T-7t,7t-7T + 7T-Bfc(T) + B fc (T)-B fc (7j)} 

< max{l - 7 T , 1 - -y T + It - B k (T)} 

< 1 - 7 T + max {0, 7T - B k (T)} , 
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hence 

sup \B k (t)- lt \<l- lT + \B k {T)- lT \ for all T>0. 
te[T,oo) 

To go further we will use the following result which resembles Dini's theorem. 



(4.26) 



Lemma 4.4 Let T > and f £ C([0, T]) &e a non- decreasing function. If (fk)k>i * s a sequence 
of non- decreasing functions on (0, P) such that f k — ► / in -^(O, T) as k — > oo, i/ien — ► / 
L°°(<J,T-5) as k -> oo /or every <5 G (0,T/2). 



Let T > rj. By (|4.20|) and Proposition ^. 31 (Bk)k>i is a sequence of non-decreasing functions that 
converges to the continuous and non-decreasing function f i— ► 7 t in L^O, T + rj) a.s. and we use 
Lemma H3] to conclude that 



lim P 



sup \B k (t)-7t\ >e/2 
te[ v ,T] 



= 0. 



(4.27) 



We now infer from (jOg} and P~2"T|) that 



lim sup P 

k — >oc 



sup \B k (t) -7t| > e 
te[i7,oo) 



sup |-Bfe(t) - 7t| > e 

tg[r(, oo) 



< P 



sup |B fc (*)-7 t | + |l-7r| >e 
teh,T] 



< P [1 - 7T > e/2] + P 

< P[l -7t > e/2] . 



sup |fl fc (t) -7 t | > e/2 
te[j?,T] 



The above inequality being valid for any T > n, we may let T — ► oo to deduce (I4.25|) since — ► 1 
as T — > oo a.s. □ 



We finally turn to the proof of the last statement of Theorem 

Proof of Point 3 of Theorem 12.51 Here (a n ) n >i is a sequence of positive real numbers such 
that a n — > oo and a n /m n — > 7 € (0, 1) as n — * 00. We consider the limit (/J.t)t>o of a subsequence 

|(^" fc:a " fc ) t >o j in the sense that a.s. 



lim sup d(fi t ' ™ fc , fi t ) — for all P > , 

fc-kOO [ 0jT ] 



(4.28) 



the existence of such a limit following by the same arguments as (|4.7[) . 
We then introduce 

Pi := inf {t > : ((i a {dx) - fi t {dx), x) > 7} , 

and aim at showing that (|Ut)te[o,Ti) solves a.s. the Flory equation (P). 
For n > 1, we set 

T? := inf{t > ; (ju™'° n , l( a „ i0o) ) > 0} , 

which represents the first time that a particle of size exceeding a n appears in the Marcus-Lushnikov 
process (/z™' a ")t>o- For each n > 1, it is clearly possible to build a classical Marcus-Lushnikov 
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process (/i™' m ") t >o (i-e. without cut-off) such that [i™' m ™ = ^"'° n for t £ [0, T™] a.s. In particular 
we have also T" = inf{< > : (/x"' m ™ (fir), l( 0n ,oo)) > 0} a - s - Denoting by M™ (t) the size of the 
largest particle at time t in the process (/z™' m ") t >o, we clearly have T™ = inf{t > : M"(t) > a n } 
a.s. 

By the tightness result of Proposition 12.41 we may assume that, up to extracting a further sub- 
sequence (not relabeled), (^™ fc ' m " fc )t>o converges to (yt)t>Q in B([0, oo), A4y ). By the Skorokhod 



(4.29) 



representation theorem and Lemma l4.2[ we may assume that a.s. 

k ,v t ) = for all T > 0. 



By Theorem 



we claim that 



lim sup d(n t 

fe^oo [ 0;T ] 



Point l-(i), (y t )t>o solves a.s. the Flory equation (F). Introducing 
S\ := inf {t > : ([io(dx) — v t (dx),x) > 7} , 



lim -T" fe | > 



for all e > . 



(4.30) 



Taking (|4.30[) for granted, we deduce that fi t — v t for t £ [0, Si) a.s. since ^™' m " = /i™' a " 
for t £ [0,T"] a.s. for all n > 1. This implies that Si — T\ a.s., because the subset {it £ 
M.^; (fj,o(dx) —ir(dx),x) > 7} = {n £ M^; (ir(dx),x) < 1 — 7} is closed in endowed 
with the vague topology, and both t 1— *■ [J4 and t 1— ► Vt are a.s. continuous for that topology by 
Lemma W% Therefore, (Mt)te[o,Ti) solves a.s. the Flory equation. 

We are left with the proof of (|4.30|) . To this end we will use ((2T5| and (|2. 10[) (with the weak 
limit (vt)t>v of the classical Marcus-Lushnikov process ( /j^ k ' m " k ) ). Introducing the (random) 

V / k>l 

gelling time S ge i of (v t )t>o given by 

S ge i := inf{t > : (v t (dx), x) < (n (dx), x)} , 

we recall that a.s. the map 1 1 — > (is t (dx),x) is constant and equal to 1 on [0, S ge i) and continuous 
and decreasing on (S ge i, 00) by Proposition ^. 31 In the proof of (|4.30p . we have to handle separately 
the events S ge i < Si and S ge i = Si, the latter being not ruled out a priori due to the possible 
discontinuity of 1 1 — > (u t (dx), x) at t — S ge i. 
Fix e > and write 

P[\Si - T" k \>e] = P[U k ] + P[V k ] + P[W k ] 



with 



w k 



{S ge i <Si< S gel + e/2 , T{ 1 " < Sx - e} , 
{S gel + e/2 < Sx , T{ lk < Si — e} , 
{S gel < Si , I** >S 1 +e} . 



First, on U k , we have a nk < M? k (Si - e) < M? k (S ge i - e/2), so that 

P[U k ] < P [M? k {S gel - e/2) >a nk ]<P 
by (I2.10p since a n /m n — ► 7 > as n — > 00. 



M^(S ge l 



/2) > Or, 
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Next, introducing r := Si — e/4 and Z := (v T (dx),x) — 1 + 7, it follows from the a.s. strict 
monotonicity of t 1 — ► (u t (dx),x) on (S ge i, 00) and the definitions of T" and Si that 



Z > a.s. and Vk C 



(1 - (v T (dx),x)) > 



~f + Z 



Let 77 > 0. For k large enough we have I7 — a nk /m nk \ < rj/2 and since r > S ge i + e/4 a.s. 



E 



> E 



sup 

te[S, ci +£/4,oo) 

M? k (r) 



M? k (t) 



m,. 



- (1 - (h> t {dx),x)) 



(1 — (v T (dx), x)) 



> E 



Z 



>E[TL Vk (Z-r)/2)] 



> 1 E [l Vk l [f , i0o) (Z)] -±E [l (0l ,)(Z)] > i P[V fc ] - 7 ? P[Z G (0,77)] 



Multiplying the above inequality by 2/77 and letting k — > 00 with the help of (|2.8[) give 

lim sup P[Vy < 2P[Z £ (0, 77)] for all 77 > . 

— >oo 

As Z > a.s., the right hand side of the above inequality converges to zero as 7/ — > 0. Consequently 
P[Vfc] — > as k -> 00. 

Similarly introducing a := Si + e and Y := 1 — (y a (dx),x) — 7, the a.s. strict monotonicity of 
t 1 — ► (v t (dx),x) on (S ge i, o°) and the definitions of T™ and Si warrant that 



y>0 a.s. and Wfe C < (1 - (v a (dx), x)) - 



Mi k (a) 



y 



Arguing as for 14, we have for all 77 > and fc large enough 



E 



sup 

t6[S flel +e,oo) 



M£(t) _ 



(1 - (v t {dx),x)) 



> ~ P[w k ]- v p[Y g (0,77)] 



We then proceed as before to deduce from (|2.8p and the a.s. positivity of Y that P[Wfe] — ► as 
k — » 00 and thus complete the proof of (|4.30p . □ 



5 Numerical illustrations 

We consider the monodisperse initial condition [Xq = Si and the multiplicative kernel K(x, y) = xy. 
Under these conditions, there is an explicit solution to the Smoluchowski equation (S) given by 



f k 



jlt(dx) :~^^c(t,k)8k{dx) with c(t, k) :- 



k>l 



fc-2 



fc! 

k k-2 



kl 



t k-l e -H 



t- x e~ k 



for t G [0, 1] . 



for t > 1 
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For the same initial condition, the Flory equation (F) has also an explicit solution given by 

h.k-2 

uJdx) := V c(t, k)S k (dx) with c(t, k) := ——t*-^-* for t > . 

fc! 

k>l 

Before proceeding to simulations, let us point out that (/i t (dx),x) = 1 for t G [0,1], while 
(/j, t (dx),x) = t*/t for t > 1, where t* G (0, 1) is the unique solution to t*e~'" — te~ l in (0, 1). Easy 
computations show that 

lnfl - t) 

Ti( 7 ) := inf{i > 0; {no(dx) - lH{dx),x) > 7} = for 7 e (0, 1) . 

7 

In Figures 1 to 5, the plain, dashed, and dotted lines represent /Lt"' a " ({2}), c(t, 2), and c(t,2), 
respectively, as functions oft. We observe that, as explained by Theorem 12. 5( 

(i) for a n *C m n , /i"' a " approximates the solution to the Smoluchowski equation, see Figure 1, 

(ii) for a n = m n , M™' " approximates the solution to the Flory equation, see Figure 2, 

(iii) for a n — r yui n with 7 G (0, 1), approximates the solution to the Flory equation until 
the time T\( 7 ), and then changes its behaviour: see Figure 3 (7 = 0.5, 7\(0.5) = 1.386), 
Figure 4 (7 = 0.8, 7\(0.8) = 2.012) and Figure 5 (7 = 0.33, T x (0.33) = 1.21). Note that 
Figure 5 shows that the behaviour of /z"' a " bifurcates at least twice on t G [0, 3]. The second 
bifurcation certainly corresponds to the time where a second giant particle with size 10 5 
appears. 
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